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Abstract — In this paper, a refinement of the weight distri- 
bution in an MDS code is computed. Concretely, the number 
of codewords with a fixed amount of nonzero bits in both 
information and redundancy parts is obtained. This refinement 
improves the theoretical approximation of the information-bit 
and -symbol error rate, in terms of the channel bit-error rate, 
in a block transmission through a discrete memoryless channel. 
Since a bounded distance reproducing encoder is assumed, the 
computation of the here-called false positive (a decoding failure 
with no information-symbol error) is provided. As a consequence, 
a new performance analysis of an MDS code is proposed. 

Index Terms — MDS code, bit-error rate (BER), block error- 
correcting code, information-bit error rate (iBER), false positive. 



I. Introduction 

A FUNDAMENTAL challenge when determining the per- 
formance of a block error-correcting code is to measure 
its bit-error rate (BER), which quantifies the reliability of the 
system. In practice, the BER estimation for a single code is 
simple, just send data and divide the errors committed among 
the total number of bits. However, it would be too costly and 
time-consuming if a comparative between several codes is 
required. Mathematical software packages for encoding and 
decoding are very limited and restricted to specific codes and 
simulations would consume a huge amount of time when 
dealing with low bit-error rates. For this reason, a theoretical 
approach to the measurement of the BER is proposed by 
several authors in the literature, see for instance |1|-|6|, |14|. 
All these papers follow this scheme of codification: let C be 
a code of length n and dimension k over the field with q 
elements, being q > 2. An n-tuple is transmitted through a 
q-ary symmetric discrete memoryless channel. In this step, 
there are two possibilities, the transmission is right or failed 
in some symbols. In a second step, the code corrects the n- 
tuple, detects an erroneous transmission but does not correct it, 
or asserts that the transmitted n-tuple is a codeword. Finally, 
there is a comparison between the encoded and decoded n- 
tuples, see Figure [T] When we run over all the possibilities in 
each step (of course, not all combinations are possible), this 
yields five disjoint cases: 
1) A correct transmission (CT), i.e., all the symbols are 
correctly received. 
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Fig. 1. Length oriented encoding-decoding scheme. 

2) A right correction (RC), i.e., some of the symbols are 
incorrectly received but the decoding algorithm corrects 
them. 

3) An error detection (ED), i.e., the number of errors ex- 
ceeds the error-correction capability of the code, the block 
is not corrected and the bad transmission is detected. 

4) A wrong correction (WC), i.e., some errors occur (be- 
yond the error capability of the code), there is a code- 
correction, but nevertheless, the encoded block differs 
from the decoded block. 

5) A false negative (FN), i.e., some symbols are incorrectly 
received but the whole block is a codeword, so, from the 
receiver's point of view, the block is correctly received. 

Cases FN and WC are called undetected errors in [6|, and it 
is proven that, for maximum-distance-separable (MDS) codes, 
the probability of an undetected error decreases monotonically 
as the channel symbol error decreases, that is, MDS codes are 
proper in the terminology of [2]. Hence, the performance of an 
MDS code is characterized by the probability that an erroneous 
transmission will remain undetected. In [1], as a performance 
criterion, the probability of an ED is added and an exhaustive 
calculus of the word- symbol-, and bit-error rate of ED, WC, 
and FN is made. 

In this paper, we propose a refinement in the calculi of 
the probability of an FN, a WC, and an ED. Consequently, 
we get a better approximation of the BER for a g-ary MDS 
code. As in the above references, we consider a bounded 
distance reproducing decoder, i.e., it reproduces the received 
word whenever there are uncorrectable errors. The underlying 
idea consists in removing the symbol errors produced in the 
redundancy part of the decoded n-tuple, that is, following the 
nomenclature of flT) , JT2J, unlike the aforementioned papers, 
we estimate the information-bit error rate (iBER), sometimes 
also called post-decoding bit-error rate. More formally, let 
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us assume, without loss of generality, that the codification 
is systematic and the first fc symbols of the n-tuples form 
an information set. Hence, following the above scheme, after 
the comparison step, if there are s errors, the symbol-error 
proportion is s/n. Nevertheless, some of these errors belong 
to the redundancy part and they will not spread in the final 
post-decoded fc-tuple. In other words, a new variable should 
be considered: the comparison between the original block and 
the final fc-tuple obtained after decoding, see Figure [2] 

Encoding 



Original fc-tuple 




n 




o 
1 








o 

3 


Target fc-tuple 




Post-decoding 

Fig. 2. Dimension oriented encoding-decoding scheme. 

Attending to this new variable, we may split the ED into 
two disjoint cases: 

3.a) A pure error detection (PED), i.e., some errors affect the 

information set. 
3.b) A false positive (FP), i.e., all errors belong to the 

redundancy part. So, from the point of view of the 

receiver, there are uncorrectable errors but, indeed, the 

post-decoded block is right. 
Then, a study of the BER should consider the probability 
of obtaining a false positive after the post-decoding process 
and, hence, the criterion for measuring the performance of the 
code should be the probability of undetected and PED errors. 

All along this paper we shall assume that the code under 
consideration is a g-ary [n, fc] MDS code in a bounded distance 
reproducing encoder, where q — 2 b for some natural b > 
and < fc < n, with minimum distance d = n — k + 1 
and error-correcting capacity up to t = \_(n— k)/2\ errors. 
Furthermore, we shall assume, without loss of generality, that 
the generator matrix is systematic in its first k coordinates. 
This is so since it is well known that, in an MDS code, any 
subset of k coordinates form an information set, see e.g. (9] 
Theorem 5.3.4]. The reorganization of these components to the 
first k does not affect our calculations and makes the text much 
more readable. The channel BER shall be denoted by p, and 
by 1 — p, the probability of a bit of being correctly received. 
The code is 2 b -ary, so the probability for a symbol of being 
correctly received is q s = (1 — p) b . Therefore, the expected 
probability of transition between any two symbols, the channel 
symbol-error rate, is the average of all possible cases, i.e., p s = 



1 ~^ 1 _i P ' > ■ Finally, if we want to know the probability for a bit 
to be erroneous inside an erroneous symbol, we simply give 
the conditional probability p b \ s = 1 _( 1 p _ p )b ■ We will denote by 
A and V the minimum and the maximum of a set of integers. 

The paper is structured as follows. In Section |TTJ we count 
the codewords of a g-ary MDS code. The number of code- 
words of a fixed weight has been calculated before (see (6| 
and [7 1); nevertheless, we need a finer approach. By means 
of arguments of Linear Algebra, we shall give a formula for 
counting the number of codewords of a fixed weight both in 
the information set and in the redundancy part. This allows us 
to calculate the iBER of an FN. In Section [Til] we shall deal 
with the decoding failures, that is, when the received word 
is corrected to a different codeword to the original one. We 
shall follow the style of 1 1 1 and calculate the words inside the 
sphere of a codeword of a given information and redundancy 
weight. We will make use of these calculations in Section [TV| 
in order to give an expression of the probability that an FP 
occurs and obtain the desired approximation of the BER of an 
MDS code. In order to make the paper self-contained, in the 
Appendix, we add some of the combinatorial formulae needed 
all along the text. 

Figures [3] to [8] have been made with Sage (lOj, using the 
Ticket # 4529: logplots.py. 

II. Counting codewords 
Let C be a g-ary [n, fc] MDS code generated by a systematic 
matrix G = (Ik\R). The aim of this section is to compute 
Aj, the number of codewords in C with i non-zero infor- 
mation symbols and j non-zero redundancy symbols, where 
i € {0, • • • , k} and j <E {0, • • • , n — fc}. We shall need the 
following lemmata. 

Lemma 1. Any square submatrix of R is regular. 
Proof: Let 

An ... ri,„_ fc ' 

R = 

V*fcl • ■ • rk, n -kj 

and let p — and 7 = {ji, ■ ■ ■ ,ji} be the sets 

of indices corresponding to the rows and columns of R that 
form an ! x I submatrix, where I < fc A (n — fc). We call 
such submatrix R pi . If det(i? P7 ) = 0, there exists a non-zero 



vector (a 1 



(«i 



, a{) € W„ such that 



■ ai) 



- 1 



(0...0). (1) 



\ r iljl r ilh ' ' ' r iljl ) 

Let a 6 be the vector in which is in the z/j-th coordinate 
for all 1 < h < I and zero otherwise. Since aG = (aIk\aR), 
there are / non-zero coordinates in the first fc coordinates of 
aG. Now, for any 1 < h < I 

r 2jh 

■ = a-indh + a 2n 23h +■■■ air njh = 
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by (TTT), so, in the last n — k coordinates of aG, there are no 
more than n—k — l non-zero coordinates. Then, aG has weight 
n — k at most. Since aG G C and the minimum distance of C 
is d = n — k + 1, we get a contradiction. Consequently, R pl 
is regular. ■ 

Given a system of linear equations, we say that a solution is 
totally non-zero if all the coordinates of the vector are other 
than zero. We say that a matrix has totally full rank if any 
submatrix has full rank. By Lemma [T] R has totally full rank. 

Let us denote by /? . the number of totally non-zero 
solutions of a homogeneous linear system, over the field ¥ q , 
with i variables and j equations whose coefficient matrix has 
totally full rank. 

Lemma 2. For any integers i,j > 1, /? . is given by the 
following recurrence: 







ifi<j 



(2) 



Proof: If i < j, there are, at least, as many equations as 
variables. Since its coefficient matrix has full rank, the zero 
vector is the only solution. Then ffj = 0. 

If i > j, since its coefficient matrix has full rank, the 
system is undetermined. Specifically, whenever we fix i — j 
coordinates, we will find a single unique solution. Then, there 
are (q — iy~ 3 solutions whose first i — j coordinates are non- 
zero. In order to calculate /? ., it is enough to subtract those 
solutions for which some of the remaining j coordinates are 
zero. For any < h < j, the solutions with exactly h coordi- 
nates being zero may be obtained choosing h coordinates from 
j, and calculating the number of totally non-zero solutions in 
a linear system of i — h variables and j equations, that is, 
(h)fi-h,j' Therefore, 



Proposition 3. For any integers i, j > 1, 



i=o ^ 7 

Proof: By Lemma [2l it is enough to check that the 
expression F? tj = £;=o (-1)* (« - l)^'"' satisfies 

(|2j. Observe that if i < j, then Ff^ — 0, since the sum runs 
over the empty set. Suppose that i > j. Since F?_ h ^ = f° r 
any h > i — j , we obtain 



J / .\ j'A(i-J-l) / .s 

Y,[ h ) F i-h,i= E { h ) F i-h,r 

h=0 v 7 h=0 v 7 



(3) 



We substitute ff_ h ,■ by i 7 ^^ j in and using ([3]), we get 

jA(*-J-l) 



i—h,j 
i—h—j—1 



E (-D 



i + 1 - 1 



j\ (j+i-i 
i 



( q - iy- h -j- 



(q - iy- h -i- 



jA(i-j-l) 

- E 

h=0 v 7 /=0 
jA{i-j-l) i-h-j-1 

= E E (-D" 

h=0 /=0 

This expression is a polynomial in (q—1). Let us now calculate 
the coefficients of (q-l) l ~ 3 ~ H for H — 0, ■ • • ,jA(i—j — l). 
We group those coefficients in which I + h = H. If H = 0, 
then i = h = 0, and therefore the coefficient of (g — l) l ~ J is 

Suppose that < H < j A (i — j — 1). The coefficient of 

(g — iy-i~ H is given by 



i. 



EC- 1 )"" 



h=0 



= EM 



(j + H-h- 1) 



7i=0 



H 



h\{j-h)\ (H-h)\{j-l)\ 



h=0 



(j + H-h-l)(j + H-h-2)---(j-h + l) 



= 0, 



where the last equality comes from Lemma 13 since (j+H 



h — X)(j + H — h — 2) • • • (j — h + 1) is a polynomial in h of 
degree H — 1. 

Because of I < i — h — j — 1, then H < i — j — 1. As a 
result the case when j<H<i — j—lis left to be analyzed. 
In such case, the coefficient of (q — l) l ~ 3 ~ H is given by 

H-hfj\f3 + B-h-l 



E(-d j 

h=a 



hj (j-l)l(H-h)l 



h=0 

t^L y (_i)^ ( i 

■(j + H-h- l)(j - if - h - 2) • • • (H - h + 1) 



0. 



where, again, the last equality follows from Lemma 13 



since 



(j + H - h-l)(j - H - h-2) ■ ■ ■ (H-h + 1) is apolynomial 
in h of degree j — 1. Then, 

jAi—j—l 



E 

h=0 



By 



F 9 = 



(9-1 



E 
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*° F h =fh- " 

We recall that A 1 - is the number of codewords of C with 

weight i in the first k coordinates (the information set) and j, 
in the remaining n — k coordinates (the redundancy part). 

Lemma 4. Let < j < n — k. If 1 < i < k, then 
/ k\ (n — k 



-4'; 



J 



E ' : 1 ( / )fi,n-k-j+r 



If i — 0, then Aj = 5ij, where 5y is the Kroneker's delta. 

Proof: The case i — is evident because the minimum 
distance of C is d = n — k + 1. Suppose that i > 1. A 
codeword with weight i in the information set and j in the 
redundancy part must be obtained as a linear combination 
of i rows of G with non zero coefficients such that exactly 
n — k — j coordinates of the redundancy part become zero, i.e. 
it is a totally non-zero solution of a homogeneous linear system 
with n — k — j equation and i variables. By Lemma [T] this 
number is given by ff n _ k _j. However some of the solutions 
counted in f q in _ k _j can turn zero some of the remaining j 
coordinates. Hence, the lemma is obtained as a consequence 
of the Inclusion-Exclusion Principle. ■ 

Remark 5. Observe that if i+j < n—k, then i < n—k—j+l, 
so Aj — 0. This agrees with the fact that the minimum weight 
of a non-zero codeword must be n — k + 1 = d. 

Let us now give a neater description of Aj. By mixing 
Proposition [3] and Lemma [4] we have 

Ai = 



k\(n-k 
V V 3 
n — k — j + 1 + h — 1 
h 



1=0 v 



-(n-k-j+l)-l 

(-!)"• 



j i— n+ fc+ j — I — 1 



l+h 



1=0 h=0 

'n — k — j + l + h — I 
h 



J 



(9-1) 



i—n-\-k-\-j — l — h 



(4) 



We proceed analogously to the proof of Proposition [3] Let 
us examine the coefficient of (q— \y- n + k +i- H m where 
I + h = H. Since h < i — n + k + j — I — 1, it is necessary 
that H < i—n + k+j — 1. We distinguish two cases. On the 
one hand, if H < j, then < H < j A (i - n + k + j - 1), 
and then, the coefficient is given by 

y^f_l\H f k \ fn-k\fj\fn-k-j + H-l 



given by 



k-j + H-l 

H - I 



H - I 



= (-i) H r 



E 



- k 



(n-k-j + H-l)\ 



- l)\ (H - l)\(n -k-j + l-l)\ 



= (-i) H r 



- k 



E 

1=0 



jl(n - k - j + H - 1) • • • (n - k)(n - k - 1)! 
l\(H - - l)l(n -k-j + l-l)\ 



= (-i) ff C 

3 

•E 



k\ (n — k\ (n — k — j + H — 1) • • • (n — k) 
V V 3 
H\ 



H{H-l)---{j + l) 

(n-k- I)! 



1=0 



i\(h - iy. (j - iy.(n -k-j + i-iy. 



/, )H f k \ ( n ~ k \ {n-k- j + H -\)---{n-k) 
1 > j J H(H-l)...(j + l) 



E 

1=0 



H\ fn-k-l 
3'1 



(-1)" 



n — k (k 



H \i 

n-k-j + H-l\fn-k + H-l 

H-l ){ 3 
t , nil — k fk\ (n — k + H — 1\ (n — k 



(-1)"' 



H H-l j \ j 

n — k\ (n — k)(n — k + H — 1)1 



H(H - l)!(n - fc)! 



=(-i) ff : 



n - k\ (n - k + H - 1)1 



={-i) h : 



j j H\(n-k-iy. 
n-k\fn-k + H-l 
3 ){ H 



(6) 



where f is given by Lemma [9] and Lemma [T0| a nd | is a 
consequence of Lemma 10 That is, by <|5j andQ), we have 
proven the main result of this section: 



1=0 



= (-!)* 
= (-1)* 



n — k 
3 

n—k 
3 



H-l 



H 

E 

1=0 

n-k+H-1 
H 



n — k — j + H 
H-l 



1 



Theorem 6. 

i— n+fc+j — 1 



H 



(5) 

where the last equality is a consequence of Lemma [9] On the 
other hand, if j < H, the coefficient of (q — iy-n+k+j-H - ls 



As a consequence of Theorem [6] we may recover the well 
known formula of the weight distribution of an MDS code as 
it appears in 16). Let {A r \ < r < n} denotes the weight 
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distribution of C. Then 



rAk 



A r — ^ A l r _ i , 



since a codeword of weight r distributes its weight between 
the first k and the last n — k coordinates. 

Proposition 7. The weight distribution of an MDS code is 
given by the following formula: 



I 1 , 



3=0 



r—j-\-l — d 



j=r-d+l 



for d<r<n,A = l and A r = Q if < r < d. 

Proof: The case r = is evident. If0<r<n — k<d, 
as we have pointed out in Remark [5 A\_ i = 0. Let us now 
suppose that r > d = n — k+1. By Theorem rol as A® = if 

j > 0, 

rAk 



£ ^r-i — 



rAk 



r— n+fc— 1 



- k 



£ (-i^l""*^" 1 ]^-!) 

H=0 



r—n-\-k — H 



We calculate ( ■) ("!•). Since r > n — k, then = 

and 



rAk 

E 



A;\ (n — k 



rAk 

E 

i=0 



/c\ (n — k 



If r > fe, then (.) = for all fe + 1 < i < r, so, for any k 
and r, 



rAk 

E 



E 

i=0 



By Lemma [9| we obtain that 

i=0 v 7 v 

so ©C^) = ("). and hence 

rAfc 



H=0 



H 



We drop (™) and develop the rest of the previous formula, 



H=Q ^ ' 

r-n+k-l , , i tj i\ 

E ( _1 ) H ( ~ 

r- n+fe-H 



E (-D 



L=0 

reordering the summands 

r— n+k— 1 r—n+k—H 

= E E 

H=0 L=0 

r — n + k — H 



L 



r-n+k-H-L 



L ,'n -k + H -I 
H 



L 



r-n+k-(H+L) 



computing the coefficient of q r n+k 3 , that is, j = H + L 

r-n+k jA{r-n+k—l) 



E E (-iy 

3=0 H=0 

r-n+k- H\ ,_„ +fc _ 3 

j-h ) q 

r—n-\-k— 1 

E (-tf- 

3=0 



n - k + H - 1 



if 



g' 



(" k + H [ ) {'' " ' '■' 11 \„ 

H=0 ^ 
r—n+k—1 



E (-i: 



r— n+k 



H=0 



n - k + H - 1 
H 



where we have split the case j = r — n + k. By Lemma 12 

for 5 = n — k, f3 = r, a = j and n = H, 



j 

E 

H=0 

and this yields 

rAk 



i? 



i=0 



r— n+/c— 1 



E (-vm vi 



3=0 

r—n+k—1 



H=0 



r-d 



E(-d 



r / ' — ' \ 7 



,r-d+l-j_ 



a \ ^ , ^ T _ d+1 (d-2 + H 



E(-!) r 



H=0 



The independent term in q of this expression can be rewritten 
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as 



r—d 



£(-*> 



r-d+l 



H=0 



d-2 + H 



H 



r—d 



d-l + H-l 
H 



( _ ir <m 

'(r -l)-(r-d)-H 
(d-l)-H 



(-1) 



r-d+l ( r 1 

d-lj' 



where the last equality is given by Lemma 12 with f) = r — 1, 
n = H, a = d — 1 and S = r — d. By Lemma [TT] with a = r, 
7 = j and k — r — d we have 



(-1) 



r-d+l 



r — d 



r-d 



4=0 



E (-D 



j— r — d-\-l 



and consequently 

rAk 



i=0 



r-d 

Em: 

4=0 



4 [ ' 1 /,»•— <*+i— 4 
J/ 



E (-D 

j— r-d+l 



This finishes the proof. ■ 
As a corollary of the above proof, we find a new description 
of the weight distribution of an MDS code. 

Corollary 8. For any r > 0, the number of weight r 
codewords of an MDS code is given by the formula 



(-1) H k + H 1 J (g_ iy-n + 1,- II 



Now, we are in the position to describe the information- 
bit and information-symbol error rate of a false negative, 
concretely, 



rjsymbol 

FN — 



rf Q n-r 



p, 



bit 



n rAk 

EE^ 

r—d i—1 

n rAk 



•Pb\s 



The iBER of a false negative can be compared for different 
codes in Figure [3] 

III. Decoding failures 

In this section, we shall make use of the calculus of the 
values A'j in order to obtain the information errors of a 
decoding failure. For simplicity, we may assume that the zero 
codeword is transmitted. Suppose that we have received an 
erroneous transmission with T\ non-zero coordinates in the 
information set and r2 in the redundancy. All along the paper, 
we shall say that the weight of this error is r\ + r 2 - If the 
received word is corrected by the code, then it is at a maximum 
distance t of a codeword. Obviously, if < n + ?*2 < t, the 




Channel BER 



Fig. 3. iBER of an FN for an MDS code over F128 of length 127 and 
dimension k. 



word is properly corrected, so we may assume that t < ri+r 2 - 
In this case, the correction is always wrong and we have 
a decoding failure. Our aim now is to count these words, 
highlighting the number of wrong information symbols and the 
errors belonging to the redundancy, i.e. the words of weight 
T\ + 7'2 that decode to a codeword c of weight ci + c 2 - 
The reasoning is as follows: for any codeword c of such 
weight, we calculate Nj- Cl ' C2 }, the number of words of weight 

fe ' (ri,r 2 )' & 

ri + 7'2 which belong to its sphere of radius t. This can be 
carried out by provoking up to t changes in c. Our reasoning 
is analogous to the one in [1|. Firstly, there is a minimum 
number of symbols that must be changed, either to zero or 
to something non-zero depending on the sign of r\ — c\ 
and r 2 — c 2 in order to obtain the correct weight. If t is 
large enough, we can use the remaining correction capacity 
to modify an additional number of symbols to zero, and the 
same number of symbols to a non-zero element of ¥ q in order 
to keep the weight unchanged. Finally, the remaining possible 
symbol modifications can be used to change some non-zero 
symbols into other non-zero symbols, without affecting the 
weight of the word. Let a = t — \c\ — r\\ — \c 2 — r 2 \ where 
\n\ denotes the absolute value of n. We may distinguish four 
cases: 



a) r\<C\ and r 2 < c 2 : In the ci non-zero information 
coordinates, ci — r\ of them should be changed to zero. 
Additionally, we also allow i% more. In the same way, on the 
c 2 non-zero coordinates of the redundancy, we must change 
c 2 — r 2 and i 2 of them. Therefore, we have the following 
number of possibilities: 



ci 

ci - r x 



C2 



C2 



> 2 



Now, we should give a non-zero value to i x coordinates 
between the k — c\ remaining information symbols, and i 2 
coordinates between the n — k — c 2 remaining redundancy 
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symbols. Thus, the number of possible changes is as follows: 



Since the changes cannot exceed t, the admissible quantities 
for i\ and i 2 satisfy 



Cl - T\ + %x + C2 — r-i + 12 + h + h < t, 



and hence 



< h + i 2 < 



Finally, we may change some of the remaining non-zero r%— i\ 
and r 2 — i 2 coordinates to another non-zero symbol. If we 
change the j\ and j 2 coordinates, respectively, we obtain 



V>- 2) TiT )U/ 2),J 



changes, where j\ and j 2 satisfy 

0< ji+j 2 <a-2i 1 -2i 2 . 
Therefore, the total number of words is the following: 



N. 



Cl \/fc-Ci 

ri — ix)\ h 



(ci,c 2 ) 
(ri,r 2 ) 

LfJ 

E 

ii=0 

LfJ 

• E 

i 2 =0 

■ £ (V 

71=0 v Ji 

a— 2i\ — 2i 2 — j 

• E 

72=0 



C2 
^2 - *2 



n — k — c 2 
i 2 



(9-1 



lH+»2 



r 2 - «2 

J2 



If we denote I = i\ + i% and J = ji + j 2 , by Lemma |9j we 
have 

N (CUC2) = 

LfJ ct-27 

= E E(9-2) J (9-1) J - 



j=o j=o 
i 



•E 

t=0 

J 

•E 

7=0 

LfJ a-27 



ci \ / c 2 \ A - cA /n - k - c 2 
n - i J \r 2 - I + i J V i J V ^ — ^ 



ri - A /r 2 - / + i 
J J\ J- J 



7=0 7=0 v 

7 , 

E 



i=0 



ci \( c 2 \ fk - ci\ fn - k - c 2 
ri-i/v r 2-^ + vV * / V ^ — ^ 



(7) 



(q 1)* 



b) r\ < c\ y r 2 > c 2 : We proceed as in the above case 
with the information symbols. So we can make 

Cl \fk-Ci 

changes in the information set. In the redundancy part, we 
must give a non-zero symbol to c 2 — r 2 + i 2 coordinates, and, 
therefore, change i 2 of C2 coordinates to zero. In that way, we 
have 

Wn-fc-c 2 \ +l2 
v «2/ \r 2 - c 2 + i 2 J 
possible changes, where 

ci - r\ + i\ + r 2 - c 2 + i 2 + i x + i 2 < t, 

and then 



< i\ + i 2 < 



Finally, changing the value of j\ and j 2 of the remaining r\— i\ 
and c 2 — i 2 non-zero coordinates, we have 



n - ii 
h 



J2 



(q - 2)*> 



changes, where 

< ji + 32 < a- 2ii - 2i 2 . 
Thus, the total number of words is 

(ci,c 2 ) _ 



Cl \/fc-Ci 

ri — i\)\ h 



(::) 




-k- 


C2 ) 






- c 2 - 


v% 2 ) 



il+V2 — C 2 +i2 



LfJ 

= E 

ii=0 

OL — 1i\ —2l2 

• E 

Jl=0 

a— 2»i — 2i 2 — ji / ■ \ 

• E ( C2 : i2 )(.-2)^. 

Again, we denote / = i\ +i 2 and J = ]\ +j 2 , and, by Lemma 
[9] we simplify the expression to 



n - ti 



AT, 



(ci,c 2 ) _ 
LfJ a-27 

£ £( g -2) J (g-l) J 

1=0 .7=0 
7 

E 

i=0 
/ 

n — i \ / c 2 — J + % 

J-J 



ci 



c 2 



+r 2 — c 2 , 

k — Ci\ [ n — k — c 2 
i 



■E 

7=0 
LfJ a-27 



EE 

7=0 .7=0 
7 

E 



r x + c 2 - I 
J 



y q -2) j ( q -i) 



I+r 2 -c 2 . 



i=0 



ci W c 2 Wfe-cA/ n-k-c 2 x 
ri - i) \ i )\r 2 -c 2 +I-i, 



(8) 
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c) 7*1 > ci y r-i < c 2 : In this case, 



N, 



(ci,c 2 ) 



( r l, r 2) 

|_§ J a-2I 



-EE ( Cl+ ,7 _/ )(?-2) J ( 9 -i) /+ -- ci 
E 



7=0 .7=0 

' Ci \ / C 2 



i=0 



k — c\ \ (n — k — C2 
' I V r*2 — I + ij \ji — ci + ij \ I — i 



(9) 



N, 



d) T\> ci y r 2 > c 2 : In this case, 



(ci,c 2 ) 



Lf J a-27 



= E E ( Cl+ j _/ )(?-2) J (g-i) /+ri - Cl - 



7=0 7=0 
7 



E 

i=0 



n — k — C2 

j / V 7 - i) V ri - ci + i) \r 2 - c 2 + I - i 



(10) 



D 



( c l, c 2) 

(ri,ra)' 



for the case ci < ?*i and C2 < 7*2, is given by 



D 



(ci,c 2 ) 
(ri,r 2 ) 



j 

•E 



Lf J a-27 



1 

( c l, c 2) 

(ri,r 2 ) J=0 7=0 
7 

Ci \ / C 2 



fc — ci\ ( n — k — C2 
. ri - 7 1 \ r 2 - I + i J V 7 / \ I — i 

4=0 X ' X ' 

(ci-ri+7)[ j 

7*1+7*2-7-1 
J - 1 



+(ri -i) 



The other cases are calculated similarly and left to the reader. 
Hence, the information-bit and -symbol error rate of a wrong 
correction is as follows, see also Figure [4] 



,symbol 



WC 



n rAk 



k n — k 



EE4-.EE ( 



jv> / — 1 p'q 



r ri — r 
sis 



r—d 2=1 



n=0 V2—0 



In order to calculate the iBER in a wrong correction, we 
follow again the ideology of (T|. Indeed, the probability of 
getting a wrong bit is different if it is due to the channel 
or due to the encoder. In the first case, it is given by the 
channel BER whilst, in the second case, it is assumed to be 



1 



9-1 



Hence, we should calculate the rate of errors in the 



information set committed by the decoder, that we shall call 
£,(ci,c 2 ) rpj^ s ca j cu j at j on j s similar to the one for numbers 

N (rt'Z)- Concretely, let us assume that ci < 7*1 and c 2 < r 2 . 
For each 7, 1, j, and J, in equation (|7]), the number of changes 
due to the decoder is c\ — 7*1 + i + j. In order to simplify the 
resulting expression, note that 



(ci - 7*1 + i) 

(ci - 7*1 + i) 
J 



7*1 + 7*2 - / 
J 

f'l + r 2 - I 
J 



Efri -i\fr 2 -I + I s 



T\ — i (v\ — i — 1\ ( r 2 — I + i 



J 



(ci - n + i) 



i-i 

?'l +7*2-1 
J 



D bit 
WC 



n rA/c 



n — fe 



i-1 n — r2— 

-i) 




Channel BER 



7-1 



(ri-i)E 



(ci - ri + 7) 



r\ — i — 1\ (t<z — I + i 

J-l-K 



r\ + 7*2 - I 
J 



+ (n - i) 



7*1 + r 2 - I - 1 
J - 1 



by Lemmas 10 and [9] where k = j — 1. Hence, number 



Fig. 4. iBER of a WC for an MDS code over F128 of length 127 and 
dimension k. 



IV. False positives 

As we pointed out in the Introduction, there exists the 
possibility of occurrence of an FR Up to our knowledge, this 
has not been treated before in the literature. In this section, 
we calculate the probability that a PED and an FP occur, 
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finishing our estimation of the iBER of an MDS code. As 
we noticed above, without loss of generality, we may suppose 
that the zero word is transmitted and we want to analyze the 
behaviour of the received word. Our purpose now is to count 
the words whose weight decomposes as + r, i.e. there are no 
non-zero information symbols, which are not corrected by the 
decoder. Obviously, if < r < t, the word shall be properly 
corrected, so we assume that t + 1 < r < n — k = d — 1. 
Two disjoint cases can take place: the error is detected but not 
corrected, producing an FP, or the error is (wrongly) corrected 
to a codeword. Since the total number of such words is given 
by 

n — k\ 



n—k 

E 

h=t+l 



h 



1 



it is enough to calculate the words corresponding to one of the 
two cases. We can make use of the calculations in Section Hill 
and give an expression of the words belonging to the second 
case. Indeed, the number of words of weight r with t + 1 < 
r < n— k that are corrected to a codeword is as follows: 

n kAs 
s—d i—1 

Hence, the number of false positives of weight r is given by 

r n — k^ 



FP(r) 



(q-iy-C(r) 



and the probability of producing a false positive is given by 
the following formula 



P FP = J2 FP (r)p r s q: 



=t+i 




Fig. 5. Probability of an FP for an MDS code over F128 of length 127 and 
dimension k. 



We may now give an estimation of the iBER of a PED. 
Indeed, when the received word has a weight greater than t, 
the error-correcting capability of the code, three disjoint cases 
can take place: an undetected error, an FP, or a PED. Hence, 




Channel BER 



Fig. 6. Probability of an FP for an MDS code over F128 of length 127 and 
dimension k (restricted range). 



for a given weight i\ + 12, the number of words producing a 
PED is given by 

PED(ix, h) = Q) (" ~ k ) (q - l)^ - FP(ia)- 



n k 



T {i,h-i) 



h—d i—1 

whenever i\ + %2 > t and zero otherwise. Therefore, 

k n — k 

ii— 1 £2 — 1 
k n — k 

p^d° 1 -EE PED^^yz+v^ 2 ■ - 



ii— 1 i2 — l 
k n — k 



p ped = E E PEDCix.ia^+^tf-* 1 - 



ii =1 i 2 = l 



k 
k 



•Pb\a 



Appendix A 
Some combinatorial identities 

For the convenience of the reader and in order to make the 
paper self-contained, we add the combinatorial identities that 
have been referenced all along the paper. 



Lemma 9 (Chu-Vandermonde Identity). 

(3 \ fct + 



p , 



5=0 



KJ \p — K 



Proof: See e.g. [8, page 8, Equation (3)]. 

Lemma 10. 

/ a\ fa — j\ 

a) \ K J \ K J V 7 
Proof: It follows directly from the definition. 



a\ ( a — k\ a fa — 1 
and — 

K I \ 7 / K \ K — 1 
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Channel BER 



Fig. 7. iBER of a PED for an MDS code over F128 of length 127 and 
dimension k. 




icr 2 10 1 



Proof: Direct consequence of the calculus of finite dif- 
ferences, see e.g. fT3] page 91]. ■ 
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Fig. 8. iBER of a PED for an MDS code over Fi 28 of length 127 and 
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Lemma 11. 

Proof: An easy induction on k using Pascal's rule. 
Lemma 12. 

/3\ _^//3-5-«\/(5 + /t-l 
a ) ^ V " ~ K / \ K 

Proof: jS] pg.7, Equation (2)] ■ 

Lemma 13. Let P{x) S ¥[x] be a polynomial whose degree 
is less that a. Then 



£(-ir- 5 ("p(<o=o 

5=0 



